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ABSTRACT 


This thesis presents a Markov model for analyzing the sensitivity of 
Crusader effectiveness to parametric changes in Crusader capabilities. Crusader 
activities are represented as finite-state, discrete-time Markov chains. A series of 
supporting models for some of the Markov model parameters provide for the 
desired sensitivity analysis on specified engineering chai acteristics of a Crusader. 
A Microsoft Excel V7.0 spreadsheet serves as the user interface, computational 
tool for the models. A Microsoft Visual Basic program manipulates the transition 
probability matrices then returns the results of the Markov model. Available 
measures of effectiveness include the time until Crusaders are killed by enemy 
counterfire and the number of fire missions executed and lost by Crusaders. A 
demonstration of the spreadsheet implementation shows the models can be used 
for the desired sensitivity analysis. While creating a model for Crusader was tiie 
motivation for this thesis, the model can be used to conduct similar analysis on 
current or competing field artillery systems. 
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THESIS DISCLAIMER 


The reader is cautioned that models and computer programs developed in this 
research may not have been exercised for all cases of interest. While every effort has been 
made, within the time available, to ensure that the models provide accurate results and the 
programs are free of computational and logic errors, they must be further validated and 
verified. The completion of these tasks is left for further research. Any application of 
these modes and programs without additional validation and verification is at the risk of 
the user. 
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EXECUTIVE SUMMARY 


The U. S. Army has identified a requirement for a new field artillery weapon 
system to support maneuver forces in battle in the next century. The Crusader and its 
companion vehicle, the Future Armored Resupply Vehicle (FARV) are under development 
to fulfill the requirement. Results of the initial Developmental Test and Evaluation 
(DT&E) and Operational Test and Evaluation (OT&E) were provided in the first Cost and 
Operational Effectiveness Analysis (COEA). The COEA identified Crusader mobility, 
range, rate of fire, and the FARV ammunition transfer rate (ATR), among others, as major 
drivers of the Crusaderfp combat effectiveness results. However, the COEA did not 
provide sensitivity analyses to show incremental changes in combat effectiveness for 
parametric changes in the Crusader capabilities. Continued DT&E/OT&E requires the use 
of transparent, low-resolution analytical modeling to conduct the desired sensitivity 
analysis on specified engineering characteristics of the CrusaderfP . 

The Markov model developed here, called the Markov Model for Crusader 
(MMC), provides for analysis of the productivity and lifetime of a Crusader pair; 
productivity is measured in the number of fire missions executed by Crusaders. The 
model incorporates the arrival of different fire mission types. Crusader availability, and 
death of a Crusader pair as a result of enemy counterfire. During its lifetime, a Crusader 
pair will wait for fire missions, execute survivability moves, conduct resupply and 
administrative activities, and execute fire missions until the pair is killed. The MMC is a 
low-resolution, analytical representation of the life and productivity of a Crusader. 

Enhancement of the MMC by the development of a series of supporting models for 
some of the MMC parameters provide for the desired sensitivity analysis. These models 
explicitly account for several engineering characteristics of the CrusaderfpPsKV: 
ammunition transfer rate, rate of fire, mobility, others. These supporting models provide 
the MMC the explicit capability to analyze the leverage a particular engineering 
characteristic has on Crusader effectiveness. Given the values for these characteristics, 
the calculations using the MMC yield measures of Crusader effectiveness in providing 
field artillery support to U. S. maneuver forces. 
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Sensitivity to the model parameters may be measured by several results of 
calculations using the MMC: 1) the expected time until a particular number of Crusaders 
is killed; 2) the expected number of fire missions executed by the Crusader pairs; 3) the 
expected number of fire missions not executed as a result of Crusaders pairs being killed 
by the enemy during firing; 4) the expected number of fire missions executed and lost by 
Crusaders at any finite time t; 5) the expected number and variance of the number of 
Crusader pairs surviving at any finite time t. The MMC provides for the analysis of a 
variety of operationally significant results. 

The Crusader Battalion Spreadsheet Implementation (CBSI) serves as the user 
interface, computational tool for the models. Microsoft Excel V 7.0 was used for the 
CBSI. The Microsoft Visual Basic programming language was used for particular 
operations within the CBSI. Provided the input parameters specified for the MMC, the 
CBSI returns the results of calculations using the MMC. A demonstration of the CBSI 
shows the Markov model can be used for the desired sensitivity analysis and, when 
compared to similar results of the COEA, provides “sensible” and “explainable” results 
with minimal computational effort. Additionally, the CBSI can be used to conduct similar 
analysis on current or competing field artillery systems. 

The Markov Model for Crusader and its supporting models provide a transparent, 
low-resolution, analytical representation of the life and productivity of a Crusader. The 
sensitivity of Crusader effectiveness to parametric changes in Crusader capabilities may 
be measured by several operationally significant results of calculations using the model. 
The user of the model simply provides the required input to the CBSI, the CBSI returns 
results with minimal computing effort. 
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I. INTRODUCTION 


A. BACKGROUND 

The field artillery is the “King of Battle” for many reasons. Field artillery fires destroyed 
the Confederate Army at Gettysburg and destroyed the will of Americans to fight at Kah Sahn, 
Vietnam. Field artillery fires not only destroy the enemy but also neutralize or suppress them so 
maneuver forces fight with great advantage. The field artillery system can fire at day and night, 
with no weather constraint. Few other weapons on the battlefield contribute as much as the field 
artillery system. [Ref 1] 

There are many reasons to modernize the field artillery system. In the Persian Gulf War, 
the speed of the Ml tanks and M2 Bradley fighting vehicles was limited to that of the slower 
M109/M992 field artillery system. The longer ranges of enemy artillery provide a significant 
counterfire advantage over the U.S. artillery system. The current system lacks the digitization and 
automation necessary to leverage technology against a numerically superior enemy force. The 
U.S. Army recognizes the need for a modernized field artillery weapon system, [Ref 2] 

The Crusader and its companion vehicle, the Future Armored Resupply Vehicle (FARV) 
are being developed for the U.S. Army to replace the current systems by the year 2006. The 
CrusaderlV will be designed with significantly increased capabilities when compared to the 
current field artillery system. The Crusader’s range will be 10-15 kilometers greater than that of 
the current system. Unlike the current system, the CrusaderlFAKV will be required to move at 
speeds equal to the Ml tank and M2 fighting vehicles. Advanced digital communications and 
information systems will increase the survivability of the Crusader/FARV over current systems 
by providing greater battlefield situational awareness. The advanced capabilities of the 
CrusaderlFAKV over existing systems will result in a modem field artillery system which will play 
a decisive role in future battles. [Ref 2] 

The Test and Evaluation Master Plan (TEMP) details the program of Developmental Test 
and Evaluation (DT&E) and Operational Test and Evaluation (OT&E) to be conducted on the 
Crusader/FARW. Results of the initial DT&E/OT&E were provided in the first Cost and 
Operational Effectiveness Analysis (COEA). The COEA identified Crusader mobility, range, rate 
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of fire, and the FARV ammunition transfer rate (ATR), among others, as major drivers of the 
Crusaderfp PiSy combat effectiveness results. Additionally, the COEA provided a comparison of 
combat results for a CrusaderFPP^RM meeting development requirements and a CrusaderfPP^SN 
with current system characteristics. However, the COEA did not provide sensitivity analyses to 
show incremental changes in combat effectiveness for parametric changes in the Crusader 
capabilities. [Ref 3] 

B. PROBLEM DEFINITION 

The updated TEMP requires the next COEA to include sensitivity analysis of combat 
effectiveness for parametric changes to CrusaderfPPiRN capabilities. The requirement calls for 
models to be constructed to perform the desired sensitivity analysis. The models must provide 
closed-form analytical results or provide for analysis on a user-interface spreadsheet. The intent 
of the requirement is to construct simple transparent model(s) to determine the leverage each 
capability has on CrusaderfpPsKM combat effectiveness. [Ref 3] 

C. PURPOSE AND OVERVIEW 

The purpose of this thesis is to present a transparent, low-resolution analytical model 
which provides for sensitivity analysis of Crusader effectiveness to parametric changes in 
Crusader capabilities. The specific objectives are to provide the following; 

1. A description of finite-state, discrete-time Markov models developed to conduct 
parametric sensitivity analysis for the Crusader. 

2. A description of the user-interface spreadsheet developed for the implementation of 
the model. 

3. A demonstration of the models’ use to conduct sensitivity analysis on the spreadsheet. 
The motivation for this study is model presentation and demonstration. 

The Markov models were developed in four phases. The Markov Model for Crusader 
(MMC) was developed first. Computational results of the MMC provide the values for 
parametric sensitivity analysis of Crusader effectiveness. The Markov Model for Crusader 
Availability (MMCA) was developed second. The resultant long-run availability of a Crusader 
pair serves as input to the MMC. A series of simple models to describe some of the input 
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parameters of the MMC and MMCA were developed third. The purpose of these models is to 
increase the fidelity in the computational results of the MMC by increasing the number of 
Crusader characteristics included explicitly in the models. The Crusader Battalion Spreadsheet 
Implementation (CBSI) of the Markov models was developed last. The CBSI provides a 
computational tool for the Markov models. 

This thesis is organized in a manner similar to the development of the models. Chapter II 
will describe the Crusader [¥and provide a summary of the operational and organizational 
concepts of their employment by the U.S. Army. Chapter III presents the formulation of the 
Markov Model for Crusader. Chapter IV presents a formulation of the Markov Model for 
Crusader Availability. Detailed numerical examples of the computations for the models are 
included in appendices to the chapters. Chapter V presents the models for some of the 
parameters of the MMC and MMCA. Chapter VI includes a description of the Crusader 
Battalion Spreadsheet Implementation (CBSI) of the models. This chapter includes a verification 
of the spreadsheet implementation and a demonstration of its use to conduct sensitivity analysis. 
Chapter VII provides conclusions, recommendations for model enhancements, and a brief 
summary. 
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n. CRUSADER AND FARV DESCRIPTION, ORGANIZATION, AND OPERATIONS 


A. CRUSADER AND FARV DESCRIPTION 

The Crusader is a 155mm, self-propelled howitzer that is being developed to provide 
indirect fires to support maneuver forces on the future battlefield for the U.S. Army. The Future 
Armored Resupply Vehicle (FARV) is the companion ammunition and fuel resupply vehicle for 
the Crusader. The Crusader and FARV “system” will provide maneuver forces a lethal, mobile, 
survivable, and sustainable indirect fire support asset. 

The Crusader will be designed to provide indirect fires for maneuver forces that are more 
timely and at greater ranges and accuracy than those furnished by the current field artillery 
weapon. The maximum range of the Crusader will increase the area that a field artillery unit can 
cover forward of its positions by up to 170% over the current system. Additionally, the Crusader 
maximum range will create a range parity with enemy artillery, thereby reducing the counterfire 
threat and making the Crusader more survivable on the battlefield than the current system. The 
Crusaded s increased rate of fire will enable two Crusaders to provide the same number of rounds 
on target simultaneously, as can a platoon of the current system. First round accuracy of the 
Crusader will be improved by an on-board Global Positioning System (GPS) which will provide 
continuous information on position location for the Crusader. The Crusader's ability to support 
maneuver forces will be limited only by an inability to provide it fuel and ammunition. [Ref 2] 

The FARV will be designed to provide ammunition and fuel to the Crusader where and 
when it needs it. The increased rate at which ammunition and fUel will be transferred fi-om the 
FARV to the Crusader will enhance the ability of the Crusader to shoot and move. The 
ammunition and fuel resupply system on the FARV will be fiilly automated which will reduce the 
personnel required to operate it by two over the current resupply vehicle. The increased payload 
of the FARV, compared to the current resupply vehicle, will provide for more rounds to be 
delivered to the Crusader. The FARV should significantly improve the Crusader battalion’s 
ammunition and fuel resupply responsiveness. [Ref 2] 

Arguably, the CrusaderIVs increased mobility will be the most significant 
improvement over the current field artillery system. Increased vehicle control and locomotion 
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capabilities will provide both systems the mobility to maneuver with forces equipped with the 
Ml A2 tanks and M2/M3 infantry fighting vehicles. Improved automotive systems will 
significantly increase the Reliability, Availability, and Maintainability (RAM) of the 
Crusader. Additionally, the improved automotive systems should reduce the 
CrusaderfPPiBN vulnerability to threat engagement by providing greater mobility. Appendix A 
provides a summary of some of the important performance characteristics of the Crusaderf^ AKV 
compared to the current M109M992-series field artillery systems. [Ref. 2] 

B. ORGANIZATION 

The CrusaderIVAKW battalion organization is expected to be similar to that of the current 
Ml 09/992-series artillery battalions. Each battalion will have a headquarters battery, service 
battery, and three identical firing batteries. Each firing battery will consist of a headquarters 
section, an ammunition supply section, and two identical firing platoons. Each platoon will have a 
headquarters element, a Platoon Operations Center (POC), four Crusaders, and four FARVs. 

[Ref 2] 

C. OPERATIONS 

The purpose of the field artillery is to provide indirect fire support to the maneuver forces. 
To accomplish this purpose, a command relationship is established between field artillery 
organizations and a maneuver force headquarters. The maneuver force commander then assigns 
each field artillery battalion a tactical mission. A field artillery tactical mission establishes a 
relationship between the battalion and a maneuver unit, or another field artillery unit. The four 
standard field artillery tactical missions are direct support, reinforcing, general support 
reinforcing, and general support. A detailed definition of these missions is included in Reference 
4. Assignment of each of these missions describes in detail the inherent responsibilities of the field 
artillery unit with respect to priority of calls for fire, zone of fire, requirements for fire support 
teams, liaison requirements, communications links, positioning authority, and fire planning. [Ref 
2 ] 

Regardless of which tactical mission is assigned to individual field artillery battalions, the 
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field artillery system must be able to perform seven basic tasks to accomplish the purpose of 
providing indirect fire support to the maneuver forces. Those seven basic tasks are coordinate 
fire support, acquire targets, deliver field artillery fires, communicate, move, maintain and 
resupply, and survive. These seven tasks are accomplished by a variety of personnel and 
equipment included in the battalion organization. A detailed description of how the 
Crusaderf¥ contributes to the accomplishment of the seven tasks and the operations in which 
they are performed is described in Reference 2. For the purpose of this thesis, a brief description 
of the deployment, resupply, and firing concepts follows. [Ref 2] 

1. Deployment Concept 

The Crusader batteries may be deployed in one of three ways: 1) by platoon; 2) by pairs; 
3) by individual sections. A Crusader platoon will occupy positions within a 1 to 2 kilometers 
(km) diameter platoon position area. The position area will be 5 to 10 km behind the forward line 
of friendly troops. Within the platoon position area, the Crusaders will typically operate in pairs. 
The Crusaders in a pair will occupy positions 500-750 meters (m) apart, but within sight of each 
other in order to provide security for one another. Within the platoon position area, the Crusader 
pairs will move (500-750 m), emplace, fire, displace to another position. At any time in position, 
the Crusaders will also be resupplied, perform maintenance, and conduct other activities. The 
Platoon Operations Centers (POCs) manage movement and operations within the platoon position 
areas. The platoon position areas will be moved periodically with the maneuver forces so that 
Crusaders remain within 5-10 km of the forward line of friendly troops. The battery commander 
manages movement of the position areas. The choice of employment technique generally depends 
on the mission of the forces, the enemy, and the terrain, among other things. [Ref 2] 

2. Resupply Concept 

The focal point for resupply operations will be the battalion Logistics Resupply Points 
(LRPs). The LRPs will be located 2-15 km behind the platoon position areas. At the LRP, each 
FARV will be reloaded with ammunition and fuel to meet the specific needs of the Crusader 
platoon that it supports. Once uploaded the FARV returns to a loiter position typically 1-2 km 
outside the platoon position area to await resupply instructions. [Ref 2] 

The counter-fire threat will dictate whether the FARVs remain in the loiter position or 
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move to their assigned platoon position area while awaiting orders to resupply Crusaders. The 
POC will dispatch a FARV to the Crusaders for resupply. When resupply is complete the FARV 
will return to its position in the platoon area, the loiter position, resupply another Crusader or 
return to the LRP to be resupplied itself [Ref 2] 

3. Firing Concept 

The firing concept for artillery consist of three components; 1) fire planning; 2) technical 
fire direction; 3) tactical fire direction. Fire planning is the continual process of selecting targets 
on the battlefield on which fires are prearranged to support particular phases of the maneuver 
commander’s battle plan. The purpose of fire planning is to optimize the use of fire support 
systems to integrate and synchronize their use at the appropriate time on the battlefield. Technical 
fire direction is the computation of firing data for the weapons systems. Primarily, the 
computations are performed with computers in the POC then transmitted to the firing system 
digitally. However, the Crusader has an on-board computer which will compute firing data as 
required to execute sensor-to-shooter missions. Technical firing computations include such 
parameters as the weather, firing element location, target location, and munitions to be fired, 
among other things. The result of the computations is the physical aimpoint settings on the 
particular firing system. Tactical fire direction is the process of analyzing fire requests to 
determine an appropriate method of attack. Tactical fire direction considers the target 
characteristics, battlefield geography, the fire control measures in effect, weapon’s effectiveness 
data, and the maneuver commander’s attack criteria. The result of tactical fire direction is a fire 
order that prescribes the method of fire, ammunition type and expenditures, units to fire, and time 
to attack. [Ref 4] 

An important product of fire planning is the Attack Guidance Matrix (AGM). The AGM 
specifies the standard fire orders and the priority in which particular target types are attacked. 

The battalion fire direction center and the POCs use the AGM to expedite formulation of the fire 
order for a particular mission. An example AGM is in Figure 1. 
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Attack 

Guidance 

Matrix 


Priority 

Target Description 

Method of Engagement 

Munitions 

Special Instructions 

1 

MRLs 

Platoon, 4 rounds 

DPICM 


2 

Artillery 

Platoon, 4 rounds 

DPICM 


3 

ADA 

Pair, 4 rounds 

HE/VT 

On-order priority 1 

4 

OP/CP 

Pair, 4 rounds 

DPICM 


5 

MVR/Tank 

Battery, 4 rounds 

DPICM 


6 

MVR/APC 

Battery, 4 rounds 

DPICM 


7 

Other 

Pair, 4 rounds 

DPICM 



Figure 1. Attack Guidance Matrix 
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III. THE MARKOV MODEL FOR CRUSADER 


This chapter describes a finite-state, discrete-time Markov Model for Crusader 
(MMC). The material in this section is taken from Reference 8. Computational results of 
the MMC provide the values for the parametric sensitivity analysis of Crusader 
effectiveness. 

A. MODEL DESCRIPTION 

The nature of Crusader operations lends itself to low-resolution analytical 
modeling by use of a finite-state, discrete-time Markov chain. Such a model assumes that 
a Crusader pair will be in any one of a finite number of states ( e.g., firing, moving, 
resupplying, etc. ) during discrete time periods. The probability distribution of the state of 
a Crusader pair at time t + 1 depends only on the state at time / and does not depend on 
the states previously occupied by the Crusader. A finite-state, discrete-time Markov chain 
model describes the stochastic nature of Crusader operations. [Ref 5] 

The MMC presented in this chapter provides for analysis of the productivity and 
lifetime of a Crusader pair; productivity is measured in the number of fire missions 
executed by Crusaders. The model incorporates the arrival of different fire mission types. 
Crusader availability, and death of a Crusader pair as a result of enemy counterfire. 

During its lifetime, a Crusader pair will wait for fire missions, execute survivability moves, 
conduct resupply and administrative activities, and execute fire missions until the pair is 
killed. 

Sensitivity to the model parameters may be measured by several results of 
calculations using the MMC; 1) the expected time until a particular number of Crusaders 
is killed; 2) the expected number of fire missions executed by the Crusader pairs; 3) the 
expected number of fire missions not executed as a result of Crusaders pairs being killed 
by the enemy during firing; 4) the expected number of fire missions executed, and the 
expected number lost by Crusader pairs at any finite time t\ 5) the expected number and 
variance of the number of Crusader pairs surviving at any finite time t. The variance, 
computed in the last result, is especially important because it can affect the number of 
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operational tests that will need to be conducted to assess the performance of the 
Crusader. It is seen that the MMC provides for a variety of significant results. 


B. MODEL FORMULATION 
1. Parameters 

The formulation of the MMC requires the specification of the following 
parameters: 

Define 


H{ty. the number of Crusader pairs alive at time t {t = 0, h, 2h, ...). 


Pm(t) 

rc 

t) 

hie) 


the probability a fire mission arrives in the time interval 
it,t + h\, 

the probability a fire mission is of type c(c=l, 
requiring c Crusader pairs to fire; 
the probability an alive Crusader pair is available for a fire 
mission in ( /, t + /i ]; 

the probability a Crusader pair engaged in a fire mission of 
type c is killed. 


A detailed description of the model parameters is included in Appendix B, Section 
A. The availability parameter, a(/f(t); /), is computed as described in Chapter IV. A 
model for the parameter, h/lc), is presented in Chapter V, Section B. 

2. Assumptions 

The following assumptions are necessary: 

• The number of Crusader pairs that are available for a fire mission 

during the interval (/,/' + /;] has a binomial distribution with Hit) trials 
and probability of success t). The current version of the MMC 

has t) constant. Future versions of this model will make this 

parameter time-dependent. 

• A fire mission is lost (rejected) if the required number of Crusader pairs 
for the particular type mission is not available during the time increment 
it, / + /?]. 

• If c Crusader pairs are engaged in a fire mission of type c, the number 
of pairs that are killed by enemy counterfire has a binomial distribution 
with c trials and probability of success /it(c). Note that this parameter 
must be specified for this model. In more comprehensive models, it 
will depend on the quantity and quality of the enemy counterfire force. 
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• Crusader pairs are considered to conduct all activities together. 

• Crusader pairs can only be killed when executing a fire mission. 

• The probability a Crusader transitions from state / to state / remains 
stationary over time. In this model, the state space is the number of 
Crusader pairs alive at time t{t = 0,h,2h,...). 

• The discrete time increment in the model, h, is 5 minutes. This time 
length is arbitrary. It is assumed a Crusader pair executes a fire 
mission and a survivability move in 5 minutes. 

• The probability that a fire mission arrives during time h in this model is 
independent of time. Future versions of this model will make this 
parameter, time-dependent. 

3. Formulation 

The number of Crusader pairs available at time t, { //(/); t = 0, h, 2h,...} is a 
finite-state, discrete-time Markov chain with stationary transition matrices of the form: 



1 

0 

0 

0 


^i(0 

QuiO 

0,2(0 - 

Q\H(0)iO 

P(t) = 

^2(0 

02,(0 

022(0 • 

02//(O)(O 



to 

to 

* Qh ( 0 ) H ( 0 ){ 0 _ 


fory>/ 

QiM = 0 (3.2) 


for 1 <j <i 


Qv{t)= 2 

prob 

prob the pJ'ob at least the number of Crusader 

<imvcs ^ , • j 1 r- « 

mission pairs needed lor the mission are 

requires at available 

least i- j 

pairs 




(i-/) 




prob i-j Crusader pairs executing the 
mission are killed 


(3.3) 


and 
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(3,4) 


HiO) 


9,(o=i-Za,(o 

j=i 


where i,j~\,2,...H(0), the number of Crusader pairs available at time 0. A three step 
procedure to compute the transition probability matrix, P(/), along with a numerical 
example, appears in Appendix C. Details of the computations using the MMC, to include 
a numerical example, can be found in Appendix D. 
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IV. THE MARKOV MODEL FOR CRUSADER AVAILABILITY 


This chapter describes the Markov Model for Crusader Availability (MMCA). The 
material in this chapter is taken from Reference 8. A result of the model, the long-run 
proportion of time a Crusader pair is available to conduct fire missions, is used in the 
MMC for the value of the parameter, a(//(0; t). Additionally, the results of the model, 
themselves, provide probabilistic information about the proportion of time a Crusader pair 
is available, conducting ammunition resupply, or under repair. 


A. PARAMETERS 
Let 


X{t) = be the state of the Crusader pair at time t {t = 0, h, 2h, ...). 
The successive states visited in the discrete-time Markov chain are: 


A = Crusader pair is available to execute a fire mission; 

B = Crusader pair is being resupplied with ammunition; 

R = Crusader pair is awaiting repair or under repair. 


The following parameters must be specified: 


Pnit) 

y 

Pb 

Pba 


Pr 


Pra 


rc 


the probability a mission arrives in a time interval of length 

the number of Crusader pairs alive; 

the probability a Crusader needs resupply in a time 

interval of length (/, / + /;]; 

the probability a Crusader pair completes resupply and 
returns to the available state in a time interval of length 
(t, t + h\, 

the probability a Crusader pair needs repair in a time 
interval of length {t, t + /»]; 

the probability a Crusader pair completes repair and 
returns to the available state in a time interval of length 
(t, t + h\, 

the probability a fire mission is of type c; requiring c 
Crusader pairs to execute; 


A detailed definition of the parameters is included in Section B of Appendix B. 
Models for some of the parameters are presented in Chapter V. These models explicitly 
account for several engineering characteristics of the CrusaderlV ARN. 
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B. ASSUMPTIONS 


The following assumptions are necessary: 

• The discrete time increment in the model, h, is 5 minutes. This time 
length is arbitrary. It is assumed a Crusader pair executes a fire 
mission and a survivability move in 5 minutes. 

• Crusader pairs are not subject to being killed in this supporting model. 
Crusaders are subject to death in the MMC. 

• Crusader pairs are equally likely to be assigned a fire mission requiring 
c Crusader pairs with probability of dy if c<y\ otherwise, the mission is 
lost. 

• The FARV is always available to resupply a Crusader pair. 

• The probability that a fire mission arrives during time h in this model is 
independent of time. Future versions of this model will make this 
parameter, pJit), time-dependent. 


C. MODEL FORMULATION 


The transition probability matrix for MMCA, P{y), is of the form: 


A 

P(y) = B 
R 


r 


Pn,Hcf\-\[Pb+Pr\ P.ilcn- Pb Pn,^rJ 
c=\ \yy c=i \yj c=i 


Pba 

Pra 




c=l 

Pba 

0 




0 

^-Pr. 


(4.1) 


where 


/ 






if-<l 


0 if->l 

y 


(4,2) 


Given the number of Crusader pairs alive isj, the probability a Crusader pair is available 
to conduct fire missions is aiy) = lim P(y)”^. Details of the computations with a 

numerical example can be found in Appendix E. 
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V. ENGmEERING PARAMETRIC MODELS OF CRUSADER MARKOV 

CHAIN PARAMETERS 


The usefulness of the MMCA and, ultimately, the MMC is measured by their 
capability to explicitly model engineering characteristics of the Crusaderfp AKW. This 
chapter presents a series of simple models for some of the parameters of the MMCA and 
MMC. The material in this chapter is taken from Reference 9. These models explicitly 
account for several engineering characteristics of the Crusaderfp AKW and provide for the 
desired sensitivity analysis using the Markov models. A detailed definition of the 
parameters and the engineering characteristics which describe them is included in Section 
C of Appendix B. Only the models are presented in this chapter, first, for parameters of 
the MMCA and second, for a parameter of the MMC. 

A. MODELS FOR THE MMCA PARAMETERS 

1. Model for the Probability a Crusader Pair Requires Ammunition 

Resupply (/?*) 

Let 

w/ = average number of rounds (projectiles) fired per mission by 

a Crusader pair; 

Aj = number of missions fired before survivability move by a 

Crusader pair; 

Lr = maximum number of rounds carried by a Crusader, 

OCR = percentage of Z,/? on-hand when ammunition resupply is 

executed. 

Let the number of missions fired between ammunition resupply equal: 

number of rounds fired 
before resupply 
average number of rounds 
fired before a survivability 
move 
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where L J represents the largest integer less than or equal to the value. Then the 
probability a Crusader pair requires resupply in a time interval of length (/, / + h] is 

Pi 1). 

2. Model for the Probability a Crusader Pair Completes Ammunition 
Resupply (/>*«) 


Let 

p = the ammunition transfer rate from the FARV to the 
Crusader (in units of Mh). 


Let the time (in units of h) to resupply a Crusader pair, Tr, equal 


= pT^ => 7 ^ = 


^ number of rounds ^ 

transferred to two Crusaders 


(2cCj^ 

ammunition transfer rate 

V J 


1 p J 


(5.2) 


Then the probability a Crusader pair completes ammunition resupply in a time 
interval of length (/, / + h] is: = min(l / 1). 

3. Model for the Probability a Crusader Pair Requires Repair (/?,) 

Let 

nif = the Mean-Time-Between-Essential-Function-Failures 
(MTBEFF) for the Crusader. 


Let the number of 5 minute time increments between failures equal /m/ x 12 (where 
there are 12, 5 minute increments in one hour and w/is given in hours). Then the 
probability a Crusader pair needs repair at the end of a 5 minute time increment is; 

= H {nif X 12). 

4. Model for the Probability a Crusader Completes Repair {pra) 

Let 

rur = the Mean-Time-to-Repair (MTTR) plus Administrative- 
and-Logistic-Downtime (ALDT) for the Crusader. 

Let the number of 5 minute time increments until repair is complete equal x 12 
(where there are 12, 5 minute increments in one hour and is given in hours). Then the 
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probability a Crusader pair completes repair at the end of a 5 minute time increment is; 

Pr. =1/KX 12). 


B, MODEL FOR THE PROBABH^ITY A CRUSADER PAIR IS KILLED AS A 
FUNCTION OF FIRING PARAMETERS 

Assume that if c Crusader pairs are engaged in a fire mission of type c, the number 
of pairs that are killed has a binomial distribution with c trials and probability of success 
h^c). 

Compute the probability of success as follows: 

Let 


Pacq 

probability Crusaders are detected by the enemy as a result 
of firing; 

e 

Crusader rate of fire in units of rounds per minute; 

rif 

average number of rounds (projectiles) fired per mission by 
a Crusader pair; 

tf 

average time to fire a mission is riflBva units of minutes; 

Ns 

number of missions fired before survivability move by a 
Crusader pair; 

Cfp = 

Crusader counterfire footprint; 

v. 

scoot velocity in units of meters per second; 

UR) = 

median time for enemy artillery to engage a Crusader after 
acquiring a Crusader in units of h minutes. The median 
time to engage is used because it is known data from 
previous Crusader analysis. [Ref 2] 


A detailed description of the parameters appears in Section C of Appendix B. 
First, the probability enemy sensors acquire a Crusader pair during firing is: 

( 

^~y~Pacqj = ^RAB ■ (5 3 ) 

This models the number of failures before the first successful acquisition of a Crusader as 
a result of firing. The exponent, 2Ns, accounts for the fact that each Crusader in a pair 
may be acquired independently during a fire mission. This model is a consequence of the 
assumption that the acquisitions of the Crusader pairs are independent and identically 
distributed Bernoulli trials with probability of success: pacq- 
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Second, the amount of time, in minutes, required for a Crusader pair to fire then 
move out of the counterfire footprint, T^s, is: 


('/K+ — (5-4) 

time to V—' 

fire mission time to 

move out 
of footprint 


Third, assume the time, in minutes, required for the enemy to engage Crusaders 
has an exponential distribution with mean 1 / pe • Compute Pe as follows; 


then 


H,={\IT,{RMn2) 


(5.5) 

(5.6) 


Therefore, the probability a Crusader pair is killed by counterfire is; 

0CRABP[SE<Tss} = [^-e-^‘^-yRAB ( 5 - 7 ) 

where Se is an exponential random variable having mean 1 / he- 

For a fire mission of type c requiring c Crusader pairs to fire, the probability of 
success in the binomial distribution for the number of Crusader pairs killed is therefore: 

= (5.8) 


C. 


NUMERICAL EXAMPLES 
1. Input 


The following is input for the examples; 


m ) 

Pmit) 

r\ 

ri 


Vs 

Cfp 

TdR) 


«/ 

e 


3 Crusader pairs; 

.3; 

. 5 ; 

• 5 ; 

8.3 meters per second; 
750 meters; 

10 minutes; 

4 rounds per mission; 
10 rounds per minute; 
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Ns 

CCr 

P 

rrif 

ntr 

Pacq 

Lr 


1 mission before survivability move; 

. 6 ; 

5 rounds per minute; 

34 hours; 

6.45 hours; 

• 7 ; 

60 rounds. 


A detailed definition of the input parameters is in Appendix B. 

2. Probability a Crusader Pair Requires Ammunition Resupply (/)*) 


N.SR, - 


.6(60)^ 

4(1) J 


3. 


4. 


5. 


6 . 


Pb=\l9 = .\\\\ 

Probability a Crusader Pair Completes Ammunition Resupply ( pba ) 


T = 


|^ 2(.6)6 0 


= 14.4 minutes 


14.4 minutes is 14.4/5=2.88 five minute time intervals 
Pba = 1/2.88 = .34722 

Probability a Crusader Pair Requires Repair (y?,) 

ntf =34 hours 
1/(34 X 12) = 00245 

Probability a Crusader Pair Completes Repair (/>„) 

Wr = 6.45 hours 
/?„= l/(6.45x 12) = .01292 

Probability a Crusader Pair is Killed as a Result of Firing 

«^=i-(i--7r=.9i 



750 

8.3(60 sec) 


= 1.906 minutes 


Pg = (1/ 10)(ln2) =.06931 minutes'' 
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The probability of success in the binomial distribution for the number of Crusader 
pairs killed is; 

/»^(c) = 9l(l-c-“®^"'^“^>) =11262. 
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VI. CRUSADER BATTALION SPREADSHEET IMPLEMENTATION 


This chapter presents the Crusader Battalion spreadsheet Implementation (CBSI). 
The spreadsheet serves as the user-interface, computational tool for the MMC and 
MMCA. The user-interface spreadsheet is described first. A simple verification of the 
spreadsheet is included second. The results of the spreadsheet are checked for reasonable 
trends. A demonstration of the use of the spreadsheet model for parametric sensitivity 
analysis of Crusader effectiveness is included third. A demonstration of the use of the 
spreadsheet model for comparative analysis with the Paladin is included last. 

A. DESCRIPTION 

The CBSI was constructed on a Microsoft Excel V7.0 spreadsheet. [Ref 10] The 
Microsoft Visual Basic programming language [Ref 11] was used for particular 
operations within the spreadsheet. The user of the spreadsheet provides the input data and 
selects a computational mode; the spreadsheet calculates and displays the results of the 
MMC and the MMCA. Instructional documentation for the CBSI software is included in 
Appendix F. 

1. Input 

The input to the spreadsheet is described in detail in Appendix B. The user must 
enter the data in appropriate unit (e.g., minutes, rounds per minute, etc.) in the 
spreadsheet. The spreadsheet then converts the input data to units of 5 minute intervals. 
For MMC finite-time results, the user must select the finite-time computational mode and 
then enter the number of battle hours for which results are desired. The input page is 
locked so that calculations do not occur until an appropriate computational mode is 
selected by the user. An example on the input page for the CBSI appears in Figure 2. The 
input data which appear in Figure 2 serve as the “base case” data for analysis in this 
chapter. The “base case” data are explained in Section B. 
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Crusader Battalion Spreadsheet Implementation 

* 12 Crusader Pairs (24 Total Crusaders) 

* Pair to Battalion Fire Missions 
User Provided Input: 

Probability Fire Mission Arrives: 


Last Updated: 6/15/96 


Probability Fire Mission is Target 

Type c Requiring c Crusader 

Pairs to Fire: 


Pair Mission 


Platoon Mission 




Battery Mission 








Battery Mission 

pm 


mm 





Battalion Mission 


Scoot Velocity = 


(750m/90s=8.3m/s) 


Counterbattery 


Footprint (meters) = 

.^ 

Crusader Payload = 



Median Time Until the Enemy Engages 
Crusaders with Counterfire (min) = 

Average Number of Rounds Fired per 
Mission - 

Rate of Fire (rounds/min) = 

Number of Missions Fired Before 
Survivability Move 

Resupply Ammunition at Percentage 
Expended (decimal) 

Ammunition Transfer Rate (rounds/min) 

Mean Time Between Essential 
Function Failures (hours) 

Mean Time to Repair plus ALDT 
(hours 



Probability Crusaders Are 
Acquired by Enemy Sensors 
as a Result of Firing 



Figure 2. CBSI Base Case Data Input Page 

2. Computations 

The computational code in the CBSI is accessible for read-only use. A password 
will be provided to the user which allows access to the code so that it can be modified for 
any reason. The “recompute” button invokes the first passage time calculations using 
formulas (D. 1) - (D. 11) of Appendix D. The “finite-time” button invokes the calculations 
using formulas (D. 12) - (D . 19) of Appendix D. 

Intermediate results for the MMC and MMCA are also accessible for examination 
in the spreadsheet model. For example, the resultant matrices of the three-step procedure 
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for computing the transition probability matrix of the MMC presented in Appendix C are 
displayed for the user. An example of the CBSI page for some of the intermediate 
computations of the “base case” input data is included in Appendix G. 

The finite-time calculations for the MMC require the most computational time in 
the CBSI. For a 24 hour battle time, a 486DX-66 computer requires 7 minutes, 45 
seconds to complete calculations and provide results. The Visual Basic code written to 
execute the finite-time calculations is displayed in Appendix H. 

3. Output 

The computational results of the MMC are displayed for the user after calculations 
are completed. An example of the CBSI results page for the “base case” input data is 
included in Figure 3. 


25 


Crusader Battalion Computational Results 


Expected Time There Are j Crusader Pairs Alive Given Start With 12 Pairs: 

(Time In Hours) 






12 

1.957216 

*6 

2.149041 



11 

1.906629 

5 

2.179115 

Expected Time Until All 

j= 10 

1.934535 

4 

2.209143 

Crusader Pairs are Destroyed 

9 

1.961758 

3 

2.284881 


26.34608 

hours 

8 

2.001763 

2 

2.416278 


1.097753 

days 

7 

2.098713 

1 

3.24701 




Expected Number of Missions Executed While 


Crusaders are Alive = 



89.77495 

Expected Number of Missions Lost While 


Crusaders are Alive = 



5.07094 

Long Run Availability Probabilities = 




Available 

Resupply 

Repair 

1 

0.884975 

0.072214 

0.042811 

2 

0.925687 

0.046655 

0.027658 

3 

0.947824 

0.032757 

0.019419 

4 

0.958943 

0.025776 

0.015281 

Crusaders 5 

0.966596 

0.020971 

0.012432 

Alive 6 

0.971719 

0.017755 

0.010526 

7 

0.97537 

0.015463 

0.009167 

8 

0.97692 

0.01449 

0.00859 

9 

0.979139 

0.013097 

0.007764 

10 

0.981038 

0.011904 

0.007057 

11 

0.982585 

0.010933 

0.006482 

12 

0.982534 

0.010965 

0.0065 

Finite Time Calculations: 


24 

1 Hours 

Expected Number of Crusader Pairs Alive 





1.645128 

Standard Error in Number of Crusader Pairs Alive 





1.921354 

Expected Total Number of Missions Executed 





76.35804 

Expected Total Number of Mission Lost 




10.00 


Figure 3. CBSI Base Case Data Output Page 


B. VERIFICATION 

The theoretical basis for a finite-state, discrete-time Markov chain model is well- 
documented and left to the interested reader for further study. The applicability of 
Markov theory for modeling Crusader operations was discussed in Section A of Chapter 
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Ill of this thesis. This section provides a brief graphical analysis to verify that the CBSI is 
a correct implementation of the MMC and MMCA. 

The input data for this analysis and the proceeding two sections are displayed on 
the CBSI input page in Figure 2 of Section A. This “base case” data were taken from 
two sources: 1) the Operational Test and Evaluation Concept for Crusader (OTEC) 

[Ref 2]; 2) the results of the Crusader Timeline Model (CTM) [Ref 6], Among other 
things, the OTEC document provides results of Crusader performance from 
OT&E/DT&E for the COEA. The CTM was developed by the Institute for Defense 
Analysis (IDA) for the Army Operational Test and Evaluation Division. The CTM 
provides information on the amount of time a Crusader spends participating in various 
operational activities. The data were taken from these sources because they represent 
operational results from previous analyses of the Crusader. The data or results may not be 
sufficient for future analysis with the MMC and MMCA. 

1. Sensitivity to the Probability a Fire Mission Arrives 

The expected time until all Crusaders are killed as a result of firing should 
decrease as the probability of receiving a fire mission in a 5 minute period increases. This 
trend reflects the fact that, as Crusaders fire more often, they are more often subject to 
detection, attack, and possible attrition. Figure 4, from results of the CBSI for the “base 
case” data, displays an intuitively appropriate relationship between Crusader death and the 
mission arrival rate. 



Figure 4. Sensitivity to Probability Mission Arrives 
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2. Crusader Pairs Alive Over Time 

The mean number of Crusader pairs alive should decrease as time increases. This 
is because the longer Crusaders are in battle, the more missions they will execute, and the 
more often they will be subject to attrition. Figure 5, from results of the CBSI for the 
“base case” data, displays an intuitively appropriate relationship between Crusaders alive, 
and elapsed time in battle. 

MEAN NUMBER OF CRUSADER PAIRS ALIVE PLUS OR 
MINUS ONE STANDARD DEVIATION 



Time in Hours 

Figure 5. Mean Crusader Pairs Alive: Sensitivity to Time 
3. Missions Executed and Lost Over Time 

The number of fire missions executed and lost should increase as time increases. 
This is because the longer Crusaders are in battle, the more missions will arrive, and the 
more often the missions will be executed or lost. Figure 6, from results of the CBSI for 
the “base case” data, displays an intuitively appropriate relationship between fire missions 
executed and lost, and time. 
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FIRE MISSIONS EXECUTED AND LOST 



Figure 6. Fire Missions Executed and Lost: Sensitivity to Time 

4. Conclusions 

There are several points worth remembering when analyzing results of the CBSI. 
First, Crusaders are killed by pairs in the CBSI. It is unlikely in battle that both Crusaders 
in a pair will always be killed as a result of enemy counterfire. Second, attrition of enemy 
artillery is not modeled in the MMC. As enemy artillery systems are killed, fewer are 
available to conduct counterfire missions which may result in less Crusader attrition. 
Lastly, the intensity of battle is not modeled in the MMC. The probability of a fire mission 
arrival should be non-stationary over time to reflect the intensity of battle. A smaller 
probability of fire mission arrival during appropriate phases of the battle will result in 
fewer missions being fired and, therefore, less Crusader attrition as a result of firing. All 
three cases may lead to greater attrition and fewer missions fired when compared to the 
results of other models. For example, simulations conducted to provide the results for the 
COEA modeled individual Crusaders, enemy attrition, and the intensity of battle and, 
therefore, probably provide results with less Crusader attrition and more fire missions 
executed in 24 hours than those results obtained from the CBSI. [Ref 7] Notwithstanding, 
the relationships between specified parameters and Crusader effectiveness are intuitively 
appealing and indicate “sensible” trends. 
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C. SENSITIVITY DEMONSTRATION 

The ultimate test of the MMC and MMCA and their spreadsheet implementation is 
whether or not they can be used to conduct sensitivity analysis on specified engineering 
parameters of the Crusader. This section provides examples to demonstrate the use of the 
CBSI for sensitivity analysis of Crusader effectiveness to changes in Crusader rate of fire, 
scoot velocity, and FARV ammunition transfer rate (ATR). A detailed description of the 
engineering characteristics is included in Appendix B. 

1. Rate of fire and Scoot Velocity 

The rate of fire is the maximum rate at which a Crusader can deliver fires. An 
increase in the rate of fire should increase Crusader survivability by reducing the time 
spent firing and thereby decrease the Crusader’s vulnerability to enemy counterfire. 

Figure 7, from results of the CBSI for the “base case” data, displays the relationship of 
rate of fire and the mean number of Crusaders alive after 24 hours of battle. 


MEAN NUMBER OF CRUSADER PAIRS ALIVE AFTER 24 
HOURS: SENSITIVITY TO RATE OF FIRE PLUS ONE 
STANDARD ERROR 


«, 5 



8 10 12 14 16 

Rate of Fire in Rounds per Minute 


Figure 7. Mean Crusader Pairs Alive: Sensitivity to Rate of Fire 

The Crusader scoot velocity is the rate at which it exits the counterfire footprint 
after executing a fire mission. The speed of the Crusader, among other things, affects the 
time to exit the footprint. The number of Crusaders alive after 24 hours should increase if 
the scoot velocity is increased because the Crusader is less susceptible to enemy 
counterfire. Figure 8, from results of the CBSI for the “base case” data, displays the 
relationship between scoot velocity and the number of Crusader pairs alive. 
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Figure 8. Mean Crusader Pairs Alive: Sensitivity to Scoot Velocity 


The number of Crusader pairs alive after 24 hours of battle appears more sensitive 
to Crusader scoot velocity than rate of fire. Improving the scoot velocity above the ORD 
requirement of 8.3 meters per second (750m/90sec) to 10 meters per second 
(750m/75sec) increases the number of Crusader pairs alive by 50%. Improving the rate of 
fire above the ORD requirement of 10 rounds per minute to 12 rounds per minute 
increases the number of Crusader pairs alive by only 12%. Additionally, the figures 
display a significant increase in the number of Crusader pairs alive as the scoot velocity 
increases, while the increased Crusader survival as a result of rate of fire increases appears 
negligible. However, the rate at which Crusader survival improves as a function of scoot 
velocity will decrease as the number of rounds fired per mission increases (4 rounds per 
mission in this analysis). Therefore, further analysis may be required to determine the 
contribution of scoot velocity to Crusader survival. 

2. FARV Ammunition Transfer Rate 

The ammunition transfer rate (ATR) is the rate at which the FARV transfers 
ammunition to the Crusader. The ORD requires an ATR of 5 rounds per minute (60 
rounds in 12 minutes). An increase in the ATR should increase the amount of time a 
Crusader pair is available to conduct fire missions. Figure 9, from results of the CBSI for 
the “base case” data, displays the relationship between ATR and Crusader availability. 
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PROPORTION OF TIME A CRUSADER PAIR IS 
AVAILABLE GIVE 12 ARE ALIVE: SENSITIVITY TO 
FARV AMMUNITION TRANSFER RATE 



Figure 9. Crusader Pairs Available: Sensitivity to 
Ammunition Transfer Rate 

While the figure displays an increase in the proportion of time a Crusader pair is 
available as a result of improving the ATR, the increase in terms of real time on the 
battlefield may be insignificant. 

D. CRUSADER vs. PALADIN DEMONSTRATION 

Certainly, the need will arise to compare the Crusader's effectiveness to current or 
competing field artillery systems. This section provides an example to demonstrate the use 
of the CBSI for analyzing the Crusader's effectiveness as compared to the current field 
artillery system; the M109A6 Paladin. 

The input data for this demonstration reflects the differences of the two systems. 
The data for the Crusader was discussed in Section B and not changed for this 
demonstration. The data for the Paladin were also taken from the OTEC and CTM. The 
Paladin engineering characteristics data (e.g., rate of fire, ATR, scoot velocity, MTBEFF, 
payload, and counterfire footprint) were inputted accordingly in the CBSI for this 
demonstration. Additionally, the probability that an arriving target is of a particular type 
was changed to reflect the fact that most enemy targets will require a Paladin platoon fire 
mission. This is because of the decreased lethality and slower rate of fire of the Paladin as 
compared to the Crusader, which will primarily participate in pair fire missions. All other 
data for the demonstration (e.g., probability a fire mission arrives, number of rounds fired 
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per mission, probability of acquisition by enemy sensors, etc.) remains the same for the 
Crusader and the Paladin. The CBSI input page for the Paladin is included in Appendix 
I. The input data for the demonstration reflects the engineering and tactical deployment 
differences of the Crusader and Paladin. 

1. Crusader vs. Paladin Pairs Alive Over Time 

The mean number of Crusader pairs alive over time should be higher than the 
mean number Paladins alive over time because Crusaders are more survivable on the 
battlefield as compared to Paladins. This superior survivability is due to, among other 
things, the Crusadefs higher scoot velocity (8.3 meters per second vs. 6.22 meters per 
second) and rate of fire (10 rounds per minute vs. 4 rounds per minute). Figure 10 
displays an intuitively appropriate comparative relationship between Crusaders and 
Paladins alive, and elapsed time in battle from results of the CBSI,. 



Figure 10. Mean Number of Crusader and Paladin Pairs 
Alive Over Time 


2. Crusader vs. Paladin Fire Missions Executed and Lost Over Time 
The expected number of fire missions executed by Crusaders should be higher than 
the expected number of fire missions executed by Paladins, while fewer fire missions 
should be lost by Crusaders than by Paladins. These effects are due to the fact that the 
mean number of Crusaders alive over time tends to be larger than the mean number of 
Paladins alive over time. Thus, more Crusaders than Paladins will tend to be available to 
execute fire missions. Figure 11, from results of the CBSI, displays an intuitively 
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appropriate comparative relationship between missions lost and executed, and elapsed 
time in battle. 


CRUSADER AND PALADIN FIRE MISSIONS EXECUTED 
AND LOST OVER TIME 



-4—Crusaders Bcecuted 
-0—Crusaders Lost 
Hi— Paladins Executed 
--B“ Paladins Lost 


Figure 11. Crusader and Paladin Missions Executed 
and Lost Over Time 


3. Analysis 

There are several observations from the Crusader and Paladin comparison worth 
highlighting. First, the mean number of fire missions executed by Crusaders is over twice 
as large as the mean number of fire missions executed by Paladins. The Crusaders^ 
higher rate of fire and greater survivability contribute to this result. Second, the mean 
number of fire missions not executed by Crusaders (10 missions lost) because they are 
unavailable (e.g., killed, under repair, executing resupply) is considerably smaller than the 
mean number of missions not executed by Paladins (52.2 missions lost). More important, 
the mean-fire-missions-executed-to-mean-fire-missions-lost ratio (mean missions executed 
/ mean missions lost) is considerably smaller for the Paladin (34.20/52.2 = .65) than for 
the Crusader (76.36/10 = 7.64). The Crusader, on the average, executes more missions 
and loses fewer than the Paladin. These results need to be verified by results from other 
models. Notwithstanding, the comparative relationships are intuitively appealing and this 
example demonstrates the usefulness of the CBSI of the MMC for such comparative 
analysis. 
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vn. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The models of this thesis lead to several conclusions. First, the Markov Model for 
Crusader (MMC) and its supporting model, the Markov Model for Crusader Availability 
(MMCA) provide a transparent, low-resolution, analytical representation of the life and 
productivity of a Crusader. Furthermore, the sensitivity of Crusader effectiveness to 
variations in the MMC parameters may be measured by operationally significant results of 
calculations using the model. 

Second, enhancement of the Markov models by the development of a series of 
simple models for some of the parameters provides for the desired sensitivity analysis of 
Crusader effectiveness with respect to specified engineering characteristics. These 
supporting models provide the user of the MMC with the explicit capability to analyze the 
leverage a particular engineering parameter has on Crusader effectiveness. 

Third, the Crusader Battalion Spreadsheet Implementation of the MMC and 
MMCA provides the computational tool to conduct the sensitivity analysis. The CBSI 
provides “sensible” and “explainable” results from the MMC and MMCA with minimal 
computational effort. 

Lastly, the MMC can easily be used to conduct similar analysis on other artillery 
systems and in particular to compare the Crusader to current or competing field artillery 
systems. Certainly, the need will arise to compare the Crusadefs effectiveness to the 
current field artillery system, the M109A6 Paladin. The inputs to the models easily 
accommodate characteristic data of the. Paladin. The Markov models and the CBSI 
provide a means to conduct alternative system effectiveness comparisons. 

B. RECOMMENDATIONS 

Time did not permit some desirable enhancements of the current Markov models 
or their spreadsheet implementation, the CBSI. The following are some enhancements 
proposed for further study. 
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One modification may be to let some of the parameters of the MMC and the 
MMCA be functions of time. For example, the fire mission arrivals may be modeled as a 
discrete-time, non-homogeneous Poisson process. This modification will provide for the 
intensity of battle to be modeled in the MMC and MMCA. 

The MMCA may be improved by enlarging the state space in the model. In its 
current form, the resultant Crusader long-run availability may be too high. While this is 
possible given the requirements for the Crusader, it more likely that, in maneuver battle, 
other Crusader activities (e.g., moving) will decrease the time Crusaders are available to 
conduct fire missions. Increasing the MMCA state space will improve the resolution of 
the MMCA by incorporating other activities of Crusaders into the model. 

The model for the probability with which a Crusader completes resupply does not 
include administrative time or time for the FARV to move to the Crusader fi-om its 
current position. These times can easily be added to the parametric model and will 
provide a more realistic representation of the resupply operations. 

The Visual Basic program in the CBSI should be improved to provide for a 
greater flexibility in the number of Crusaders and mission types available to the user. For 
instance, the need may arise to analyze a Division Artillery ( 3 battalions = 72 Crusaders). 
The MMC provides for analysis of any number of Crusaders and mission types. In its 
current form, the CBSI does not allow this flexibility. A maximum of 12 Crusader pairs 
(battalion) are modeled in the CBSI. This improvement will eradicate the current 
requirement that the 24 Crusaders in a battalion conduct all activities together. 

The results of Chapter VI demonstrated that the Markov models’ behaviors agree 
with intuition. However, these results cannot be considered to be a complete validation of 
the models. The results of the MMC need to be further compared to results of other 
models, and to field exercise data, if available. Additionally, the distributional assumptions 
for the parameter models need to be tested against data for the activities the represent. 
While every effort has been made, within the time available, to ensure that the models 
provide accurate results, further validation efforts are necessary to ensure complete 
confidence in their predicted outcomes. 
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C. SUMMARY 

The U. S. Army has identified a requirement for a new field artillery weapon 
system to support maneuver forces in battle in the next century. The Crusader and its 
companion vehicle, the Future Armored Resupply Vehicle (FARV) are under development 
to fulfill the requirement. Continued DT&E/OT&E requires the use of transparent, low- 
resolution analytical modeling to conduct sensitivity analysis on specified engineering 
characteristics of the CrusaderlV ASCsJ. 

The purpose of this thesis has been three-fold: 1) to describe the finite-state, 
discrete-time Markov models developed to conduct the desired sensitivity analysis on 
specified engineering characteristics of a CrusaderFFAKW, 2) to describe the user- 
interface spreadsheet developed for the implementation of the Markov models; 3) to 
demonstrate the models’ use to conduct the desired sensitivity analysis using results from 
the spreadsheet implementation of the model. 

The nature of Crusader operations lends itself to modeling by a Markov chain, and 
the Markov models developed in this thesis were demonstrated to be useful. However, 
the results are only as good as the assumptions used to create the models and the accuracy 
of the values for estimated parameters. Implementing the recommended enhancements to 
the models should increase usefulness and validity in the results obtained from the current 
models. Empirical studies should be conducted to estimate the models’ parameters. After 
doing so, the MMC and its supporting models will provide flexible, low-resolution, 
analytical models for the parametric sensitivity analysis of Crusader effectiveness. 
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APPENDIX A. COMPARISON OF PALADINfYAASW AND CRUSADERI¥AR\ 
OPERATIONAL REQUIREMENT CHARACTERISTICS 


This appendix presents some of the critical operational requirements of the 
CrusaderPEAKV and the current field artillery system, the PaladinfEAASW. These data 
were taken from the Operational Requirements Document for the Crusaderf¥ ARN and 
summarized in the OTEC [Ref 2] . Reference 2 includes a complete listing of similar data 
on the systems. Table 1 and 2 display some characteristic data for the systems. 


Critical Characteristic 

Crusader 

Paladin 

Maximum Range 

40-50 km 

30 km 

Rate of Fire (burst) 

10 rounds/minute 

4 rounds/minute 

Payload 

60 rounds 

39 rounds 

Cross-Country Mobility 

39-48 km per hour 

27 km per hour 

Scoot Velocity 

8.3 meters/second 
(750m/90s) 

6.22 meters/second 
(560m/90s) 

MTBEFF 

34 hours 

21 hours 

MTTR 

2 hours 

2 hours 

Crew Size 

3 soldiers 

4 soldiers 


Table 1. Operational Requirement Characteristics of the Crusader and Paladin 


Critical Characteristic 

FARV 

FAASV 

Ammunition Transfer 
Rate 

60 rounds in 12 minutes 

60 rounds in 35 minutes 

Payload 

130 rounds 

90 rounds 

Cross-Country Mobility 

39-48 km per hour 

31 km per hour 

Crew Size 

3 

5 


Table 2. Operational Requirement Characteristics of the FARV and FAASV 
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APPENDIX B. INPUT PARAMETERS DESCRIPTIONS 


This appendix provides a description of the parameters used in the MMC, MMCA, 
and CBSI. 


A. MARKOV MODEL FOR CRUSADER PARAMETERS 

Parameters 

Pm{t) This parameter represents the stationary probability that an enemy 

target will arrive during a 5 minute time increment (/, t + h\. The 
greater the probability, the greater the number of fire missions 
assigned to Crusaders. This probability is effected by the intensity 
ofbattle. 

rc This parameter represents the probability that a target arrival is of 

type c. A type c target requires c Crusader pairs to execute the fire 
mission. The arrival of a type c target is effected by the phase of the 
battle and the force composition of the enemy. The number of 
Crusader elements required to attack the target is dependent on the 
Crusader rate of fire and lethality, and the size, disposition, and 
importance of the target. The Attack Guidance Matrix described in 
Chapter II of this thesis prescribes the pair to battalion size 
Crusader element to engage a type c target. 

a(H(ty, t) This function represents the probability that an alive Crusader pair 
is available to execute a fire mission. In battle, this parameter is a 
function of many factors (e.g., MTEFF, resupply time, movement 
time, etc.). In the MMC, the availability is a function of the 
probability of requiring repair or requiring ammunition resupply. A 
result of the MMCA, the long run proportion of time a Crusader 
pair is available to execute a fire mission, is used as the value of this 
parameter. A description of the MMCA is included in Chapter IV 
of this thesis. 

hdp) This parameter represents the probability a Crusader pair engaged 

in a fire mission of target type c is killed by counterfire. Fire 
missions requiring more than one Crusader pair may increase the 
probability that a single pair is killed by counterfire. In battle, this 
parameter is a function of many factors (scoot velocity, rate of fire, 
enemy sensor and counterfire capability, etc ). 
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B. MARKOV MODEL FOR CRUSADER AVAILABILITY PARAMETERS 


Parameters 

Pmif) 


y 


Pb 


Pba 


Pr 


Pra 


This parameter represents the stationary probability that an enemy 
target will arrive during a 5 minute time increment (r, t + h\. The 
greater the probability, the greater the number of fire missions 
assigned to Crusaders. This probability is effected by the intensity 
of battle. 

This parameter represents the number of Crusader pairs alive This 
parameter is exogenous to the MMCA model since there is no 
attrition. 

This parameter represents the probability a Crusader pair needs 
ammunition resupply at the end of a 5 minute time increment. The 
larger the probability, the more time a Crusader pair spends being 
resupplied and, consequently, the less time a Crusader pair is 
available to execute fire missions. 

This parameter represents the probability a Crusader pair completes 
ammunition resupply at the end of a 5 minute time increment. The 
larger the probability, the less time a Crusader pair spends being 
resupplied and, consequently, the more time a Crusader pair is 
available to execute fire missions. 

This parameter represents the probability a Crusader pair needs 
repair due to an essential function failure at the end of a 5 minute 
time increment. The larger the probability, the less time a Crusader 
pair spends being repaired and, consequently, the more time a 
Crusader pair is available to execute fire missions. 

This parameter represents the probability a Crusader pair completes 
repair at the end of a 5 minute time increment. The larger the 
probability, the less time a Crusader pair is under repair and, 
consequently, the more time a Crusader pair is available to execute 
fire missions. 


This parameter represents the probability that a target arrival is of 
type c. A type c target requires c Crusader pairs to execute the fire 
mission. A detailed description of this parameter is included in 
Section A of this appendix. 
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C. INPUT PARAMETERS 

Parameter 

This parameter represents the average number of rounds 
(projectiles) fired by a Crusader on a fire mission. The 
number of rounds fired is determined by the desired effects on the 
target and the lethality of a particular type round. More rounds are 
required to destroy a target then to suppress it. Standard data 
tables provide information on the effects of particular munitions on 
particular target types and the number required to attain a desired 
effect with a given probability. The ammunition available may 
limit the number of rounds fired. Ultimately, the guidance on the 
number of rounds to fire is published in the Attack Guidance 
Matrix before an operation. 

0 This parameter represents the maximum rate at which Crusaders 

fire rounds in units of rounds per minute. The greater the rate of 
fire, the more rounds that can be fired in a specified time. This rate 
can be sustained for 5 minutes at which point it decreases because 
of the temperature of the Crusader gun. The rate is an engineering 
parameter constrained by the manufacturing characteristics of the 
gun. 

Ns This parameter represents the number of independent fire missions 

a Crusader executes before it conducts a survivability move out of 
the enemy artillery counterfire footprint. The more missions 
executed before movement the more likely a Crusader is subject to 
detection and attack. The artillery commander determines this 
parameter based on weighing the need to provide continuous fires 
(uninterrupted by movement), the enemy detection capabilities and 
counterfire threat. 

Lr This parameter represents the maximum number of rounds carried 

by a Crusader. It is a physical manufacturing characteristic of the 
Crusader. 

(tR This parameter represents the percentage of the maximum number 

of rounds carried by a Crusader (Lr) at which time resupply is 
executed by a FARV. The percentage is determined by the artillery 
commander and is effected by the battle plan and the amount of 
available ammunition. 

p This parameter represents the rate at which the maximum number 

of rounds carried by a Crusader (Lr) can be resupplied by a FARV 
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in units of rounds per minute. The rate is determined by the 
manufacturing characteristics of the FARV. 

This parameter represents the mean-time-between-essential- 
function-failures (MTBEFF) in units of hours. An essential 
function failure prevents a Crusader from performing its mission 
(e.g., cannot deliver fires, cannot move, etc.). The MTBEFF is 
effected by strictly manufacturing characteristics. 

This parameter represents the mean-time-to-repair (MTTR) a 
Crusader plus the administrative-and-logistics-downtime (ALDT) 
in units of hours. The MTTR is the time required to repair an 
essential function failure. The ALDT is the time a Crusader spends 
waiting for repair personnel or required repair parts. 

This parameter represents the counterfire footprint; the distance (in 
meters) a Crusader must move (in any direction) after firing to 
decrease to an appropriate level the probability of attrition by 
enemy counterfire. The distance is based on the lethality and 
accuracy capabilities of enemy artillery and sensors. The smaller 
the distance, the less likely a Crusader is attrited by enemy 
counterfire. 

This parameter represents the scoot velocity; the rate at which a 
Crusader exits the counterfire footprint after conducting a fire 
mission in units of meters per second. The greater the scoot 
velocity, the less likely a Crusader is subject to attrition. The scoot 
velocity is dependent on the size of the counterfire footprint and the 
time required to move out of the footprint. The scoot velocity is 
computed as follows: 

Cfp / time to exit footprint (D 1) 

The scoot velocity is affected, explicitly, by the speed and 
displacement time of the Crusader. 

This parameter represents the probability that, as a result of firing, a 
Crusader is detected by enemy sensors. This probability is effected 
by the firing signature and rate of fire of the Crusader, and the 
quantity and quality of enemy sensors. 



Te(R) This parameter represents the median time for enemy artillery to 

engage a Crusader after acquiring the Crusader. This time begins 
to accumulate after enemy sensors detect a firing Crusader element 
The smaller the time, the more likely enemy counterfire will result 
in Crusader attrition. The time is dependent on the quality of the 
enemy fire support system. 
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APPENDIX C. THE TRANSITION PROBABILITY MATRIX FOR THE 
MARKOV MODEL FOR CRUSADER 

This appendix provides a three-step procedure to create the transition probability 
matrix, P{f)^ for the Markov Model for Crusader. A numerical example is included to 
demonstrate the procedure. 


A. THE MATRIX MKILL 

Compute the square matrix; 

= P{j Crusader pairs are killed|mission is of type c} = (c)^ [l - 

j = 0,1,-c andM^(cj) = 0 forj>c. (C.l) 

The number of rows and columns are the number of different mission types. Each entry of 
the matrix is the probability j Crusader pairs are killed of the c that are assigned to a 
mission (binomial distribution). The probability of success, h^c), is described in Chapter 
V. 




M,iC-\) 


M,(2-2) 

M,{C-2) 


M,{2-C) 

M,(C,C) 


(C.2) 


B. THE MATRIX KKILL 

Compute the matrix Kyit) where 

Kij(t) = Pjy Crusader pairs are killed]/ Crusader pairs are alive |. 

Ify>/ 


and if 1 <7 < / 


K,j(t) = 0 


(C.3) 


mm(C,0 I 

^(0= i: (c.4) 

c-j k~c 

The number of rows of the matrix is the number of Crusader pairs alive at i = 0. The 
number of columns of the matrix is the number of mission types, c. 
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Additionally, compute: 

Bq 0;0 = P\zQTo crusader pairs are killed]/ are available and a mission arrivesj 

= 1- (C5) 

;=i 


C. THE TRANSITION PROBABILITY MATRIX 

Compute the entries for the transition probability matrix, P(t), as follows; 



0,(0 = 



(C.6) 

! 

II 

cS 

m(0)+p™(0^o0;0 

(C.7) 


qAt)-- 

= i-ZQs(0. 

7=1 


(C.8) 

1 

0 

0 

0 


(0 

0n(O 

0.2(0 - 

Q\H{{3) (0 


^2(0 

02,(0 

022(0 ••• 

QlHiO) (0 

(C.9) 


Qh{0)\ (0 

0//(O)2(O ••• 




D. NUMERICAL EXAMPLE 

1. Input 

Pm{t) = . 3 ; 

r\ = .5 ; 

r2 = .5 ; 

a{H{t) = 1; /) = .928995 computed in Appendix E; 

a{H{t) = 2;t) = .897144 computed in Appendix E; 

a(H{t) = 3; f) = .928995 computed in Appendix E; 

hiXc) = . 11262 computed in Chapter V; 

H(0) = 3 Crusader pairs at time t = 0. 

A detailed definition of the input parameters is included in Appendix B. 

2. Computations for the Matrix MKill 

^;^0;l) = (l)ll262’[l-.11262f''^ = .11262 
M^(2;l) = (^).11262'[l-.11262]^"‘'^ = .19987 
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M^(l;2)=: 0 

M^.( 2 ; 2 ) = (^).11262^[l-.11262]^^‘^^ = .01268 

, '.11262 0 1 
[.19987 .01268 

3. Computations for the Matrix Kkill 

K,,(t)= .5[i-([).928995°(1-.9289951)''“|m^(1,1)= .052312 
K,,(t)=0 

^2.(0 = .5[l -(^).897144°(1-.897144 )^"“]m,,(1,1) 

+.5[l - ((I ).897144° (1-.897144)'-° + ([ ).897144' (1-.897144)'-' 
*M^(2,\)= .580434 

^ 22 ( 0 = .5[i-(^).897144°(1-.8971441)^-°]m^(1,2) 

+.5[l-((^).897144“(l-.897144)"'“+(^).897144'(l-.897144)^'' 
*M^(2,2)= .005103 
(0 = -Sjl - (^).928995®(1-.928995 ])'-‘’]m^. ( 1 , 1 ) 

+.5[((^).928995^(1-.928995)'-' + (^).928995^(1-.928995)"-")] 
*M^(2,l)= .598495 

^32 (0 = - (o)-928995“(1-.928995)'-°]m^ ( 1 , 2 ) 

+.5[(^).928995'(l-.928995)'-' + (^).928995'(l-.928995)'-'' 

* Ml, (2,2)= .006249 

'.052312 0 

K..(t)= .580434 .005103 
.598495 .006249 

4. Computations for the Transition Probability Matrix 

000 (0= 1 

a, (/) = (l-.3)+.3[l - K, , (/)] = .98431 
022 (0 = (l-.3)+.3[l - (^21 (0 + ^22 (0)] =.81858 
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Q„(l) = (l-.3)+.3[l (0 + X„(0)] = .81858 

0^^(t)=3K^, = .17413 
023(0= 0 

.00187 

Q,,it)=3K,,= .17955 

0.0 (0 = 1 - (0,, (0 + 0,2 (0 + 0,3 (0) = 01569 

020 (0 = 1- ( 02 , (0 + 022 (0 + 023 (0) = 00729 

030 (0 = 1- ( 03 , (0 + 032 (0 + 033 (0) = 0 

1 0 0 0 

.01569 .98431 0 0 

.00729 .17413 .81858 0 

0 .00187 .17955 .81858 
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APPENDIX D. COMPUTATIONAL FORMULAS FOR THE MARKOV MODEL 

FOR CRUSADER 


This appendix presents the computational formulas of the MMC described in 
Chapter III. The computational formulas in this section were taken from Reference 8. 
These formulas provide for parametric sensitivity analysis of Crusader effectiveness. A 
numerical example is included to demonstrate the computations. 

A. COMPUTATIONAL FORMULAS 

1 . The Expected Time Until j Crusader Pairs are Killed 

Let 

Na— the number of time intervals during which there are a Crusader pairs 
alive. 

Then if { H{ty, t = 0, /2, 2/;,. . . } is a Markov chain having stationary transition 
probabilities, 


£[AfJi/(O) = /] = (/-0-; (D.l) 


where / is an identity matrix with the same dimensions as Q and Qy is defined in (3.1). 
Notice that Qy is the submatrix of the transition probability matrix for the MMC. 

Let the time until at least j Crusader pairs are killed be 
Tk (y) = min{ t > 0; //(O < H{(S) - j) . 

Then the expected time there are (j + 1) or more Crusaders alive given /, the initial 
number of Crusader pairs, is: 


E[T^ (J)\m) ='] = £[a^, + N,., + l^?(0) = /] 

= 

k= ;■+! 


(D.2) 
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2. The Expected Number of Missions Executed while Crusaders are Alive for a 
Markov Chain Model with Stationary Transition Probabilities 


Let 

S(t) = the number of fire missions executed during (0, t\ 

R(t) = the number of fire missions lost (rejected) during (0, r] because not enough 
Crusaders were immediately available to fire the mission. 


Then 


(a; t+h) = + h)\H{Q) = a] = 


P.t 

prob 
mission 
arrives 
in (0,h] 

V._ 


3 

prob _ _ i 

rnission ^ Crusader 

pairs are available 


expected number of missions 
executed during (0, h] 


+Z 

^ ~ ^ prob j expected number 

pairs are of missions 

alive after executed during 

the mission the remaining 

time intervals 
given j pairs are 
alive 


(D.3) 


where if c > a 


j^[:)cciay(l-a{a)r-‘ = 0. (D.4) 


Let 


(a; oo) - lim E[Sit)\ //(O) = a] 


= Pm 






;=i 

Let /Wj(qo) be a column vector withy* entry ms(j; oo). Then oo) can be rewritten as 

(;;«>) = /.(«)• (D.6) 
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Let / be an ( H(0) x H(0) ) identity matrix. Then ( 1-Q is a square matrix with H(0) 
rows. Let/j be a column vector with a entry fs(d). Then ( 1-Q Yfs is a column vector with 
/th entry 

(d.7) 


Therefore, 

- 2 ]'/.](«); (D.8) 

the entry of the column vector is the expected total number of fire missions 

executed given H(Q) = a. 

3. The Expected Number of Missions Lost while Crusaders are Alive for 
a Markov Chain Model with Stationary Transition Probabilities 

Let 


R(a) = the number of missions rejected while there are a Crusader pairs 
alive; 

N(a) = the number of time intervals there are a Crusader pairs alive. 


Then 


£[;?(a)|A^(a)]= A^(a) 


P.jl n, Z(,”)a(<i)'(l-a(a))"- 


prob prob 

mission mission 

arrives of type c 

in (0,h] 


prob fewer than c Crusader 
pairs are available 


(D.9) 


Thus, 


C c-\ 


£[/?(«)] = E[E[i?(a)|Ar(a)]]= )«(«)'(I - «(«))'’“'. (D.IO) 

c=\ 1=0 

Therefore, the total expected number of missions lost while there are Crusaders alive is: 


ff(0) 




^J=l 


H( 0 ) 


C c-1 


a=] c-\ 1=0 
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4. Expected Number and Variance of Crusader Pairs Alive at Time t 
The n-step transition probability matrix for (3.1) of the MMC is: 

= P{H{t) = j\H{Q) = /} = [/>(!) X P(2)x...xP{l)\. (D. 12) 

Then, the expected number of Crusader pairs alive at time t is: 

£[i/(()|ff(0) = (D 13) 

J 

In order to determine the variance, first compute: 

E[H{tf\H{0) = /] = (D 14) 

j 

Then the variance in the expected number of Crusader pairs alive at time t, by definition, 
is: 

Var[H{,t)\H{0) = /] = E[H{tY\H{Q) = /]- £[//(/)|i/(0) = /f . (D. 15) 

5, The Expected Total Number of Fire Missions Executed and Lost 
During ( 0, t ] for Markov Chain Model with Non-Stationary 
Transition Matrices 

First, the probability a fire mission is executed or lost during {t, / + /?] must be 
computed as follows: 

J '' , . _ 'prob prob k-c 

prob J Crusader pairs mission mission^ , . , Z 

are alive at time t arriver^ of type ^ Crusader pairs 

are available 
c 

= Z P{n(t) = J\m = (D. 16) 

J 

and 
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PlosM^t) = 

J ' -—^-' nr ^ <^1 *=0 


prob j Crusader pairs 
are alive at time t 


prob 

mission 

arrives 


prob 
mission 
of type 


prob fewer than c Crusader pairs 
are available 


^ - 7l^(0) = i}Piosht) ■ (D17) 


Then the expected total number of missions executed during (0, / ] given //(O) = / is 

t-h 

fnexec 0; 0 = Z Pe^c s) (D. 1 8) 

and the expected total number of missions lost during ( 0, / ] given H(0) = i is 

t-h 

f^hsi O'; 0 = Pios, (j ; -y) • (D • 19) 

s=0 


B. 


NUMERICAL EXAMPLE 
1. Input 


Pnit) 

n 

ri 

a(//(0=l;0 - 
a{H{t) = 2-t) = 
a(//(0 = 3;0 = 

Me) 

m 


•3; 

.5; 

.5; 

.928995 computed in Appendix E; 
.897144 computed in Appendix E; 
.928995 computed in Appendix E; 
.11262; 

3 Crusader pairs at time t = Q. 


2 . 


Expected Time Until J Crusaders are Killed Given Start with Three 
Crusader Pairs 


{i-QT 


63.734 0 0 

61.1738 5.51207 0 

61.2002 5.455 5.512 


£[4(1)|M0) = 3] = 5.512 


£[r^(2)|M0) = 3] = 5.512 + 5.455 = 11.0575 


£[7’^(3)|M0) = 3] = 61.2002 +11.0575 = 72.2577 


55 


The expected time until all three Crusader pairs are killed is: 72.2577 five minute time 
increments (approximately 6 hours). 

3. Expected Number of Missions Executed 

/,(1) =.3[.5(l - (j,).928995'’(]-.928995)''°)] = .139349 

/, (2) =.3[.5(l - (0 ).897144 “ (1- .897144))+.5(' ).897144' (1-.897144)] = .269143 
/X3)=.3* 

[. 5(1 - (^ ).928995'’ (1-.928995)''® )+.5(( ^ ).928995' (l-.928995)'‘' + (^ ).928995' (1-.928995'- 
= .297785 

'.1393491 [8.881405 
(/-0'7.(7) = (/-0'' -269143 = 10.00806 

.297785] [ 1 1.63786 

The expected number of missions executed = 11.63786. 

4. Expected Number of Missions Lost 

£[7?(l)|/7(0) = 3 ] = [/ -0]‘'.3[.5(j,).928995°(l-.928995)'"'’+.5(l)] = 9.831858 . 

£[R(2)|//(0) = 3] = [/-2E 

*.3[.5(^).897144‘’(1-.897144)'‘‘’+.5((^).897144“(1-.897144)"''’ +(f).897144'(1-897144)""')] 

= 168332 

£[/?(3)|//(0) = 3] = [/-Q];; 

*.3[i(^).928995‘’ (l-.928995)"-“ +.5((^).928995'' (1-.928995)"-" +(J)928995' (1-.928995))] 

= 01221 

" 3 " 

E 2^7?(a)|/7(0) = 3 = 9.831858+.168332+.01221 = 10.012 

_a~] _ 

The expected number of missions lost is 10.012. 
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5. Expected Number of Crusaders Alive After 20 Minutes of Battle (4 
Time Intervals of 5 Minutes in Length) 

10 0 0 
.01569 .98431 0 0 

.00729 .17413 .81858 0 

0 .00187 .17955 .81858 

10 0 0 
.031134 .968866 0 0 

.01599 .313937 .670073 0 

.001338 .034636 .293952 .670073 

1 0 0.0 
.046335 .953665 0 0 

.0258 .425691 .548509 0 

.004025 .086532 .360935 .548509 

1 0 0 0 

.061298 .938702 0 0 

.036748 .514524 .448998 0 

.008014 .14905 .393939 .448998 

£[//(20 min )|/7(0) = 3 ] = 0(.008014) +1(.14905) + 2(393939) + 3(448998) = 2.28392 
The expected number of Crusader pairs alive after 20 minutes is: 2.28392. 
Far|//(20min)|//(0) = 3j = [o^ (.008014) + F (.14905) + 2^ (.393939) + 3" (.448998)] 

-2.28392^ =.54950 

The standard error in the expected number of Crusader pairs alive after 20 minutes of 
battle is: -\/.54950 =.74128 Crusader'pdirs. 

6. The Expected Total Number of Missions Lost in 20 Minutes 
For this section and the proceeding section in this appendix the cumulative 

binomial distribution computational form will be represented by the following notation: 
bin(number of successes, number of trials, probability of success). Then the expected total 
number of fire mission lost during 20 minutes of battle is: 
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p^(0) = (.3)[.5bin(0,3,.928995)+.5bin(l,3, .928995)] =.002215 

/?^(5min)=.0(.3)[.5(l)+.5(l)]+.00187(.3)[(.5)bin(0,l,.928995)+.5(l)] 

+.17955(.3)[.5bin(0,2,.897144)+.5bin(l,2,.897144] 

+.81858(.3)[.5bin(0,3, .928995)+.5bin(l,3, .928995)] = 00765 

/?^(10min) =.001338(.3)[.5(l)+.5(l)]+.034636(.3)[(.5)bin(0,l,.928995)+.5(l)] 
+.293952(.3)[.5bin(0,2,.897144)+.5bin(l,2,.897144] 

+.670073(.3)[.5bin(0,3, .928995)+.5bin(l,3, .928995)] =.02013 

(15 min) =.004025(.3)[.5(l)+.5(l)]+.086532(.3)[(.5)bin(0,l, .928995)+.5(l)] 
+.36093 5(.3)[.5bin(0,2, .897144)+.5bin(l,2, .897144] 

+.548509(.3)[.5bin(0,3, .928995)+.5bin(l,3, .928995)] =.02746 

(20 min) = (0) + (5) + p^ (10) + p^ (15) =.07739. 

The expected total number of missions lost in 20 minutes is .07739. 

7. The Expected Total Number of Missions Executed in 20 Minutes 

PsiG) = (.3)[.5(1 - bin(0,3,.928995))+.5(l - bin(l,3,.928995))] = 297785 

^^(5) = 0(.3)[.5(0)+.5(0)]+.00187(.3)[(.5)(l - bin(0,l,.928995))+.5(0))] 
+.17955(.3)[.5(1 - bin(0,2,.897144))+.5(l - bin(l,2,.897144))] 
+.81858(.3)[.5(1 - bin(0,3,.928995))+.5(l - bin(l,3, .928995))] =.29235 

^^(lO) =.001338(.3)[.5(0)+.5(0)]+.034636(.3)[(.5)(l - bin(0,l, .928995))+.5(0))] 
+.293952(.3)[.5(1 - bin(0,2,.897144))+.5(l - bin(l,2,.897144))] 
+.670073(.3)[.5(1 - bin(0,3, .928995))+.5(l - bin(l,3, .928995))] = 28348 

Ps(\S) =.004025(.3)[.5(0)+.5(0)]+.086532(.3)[(.5)(l - bin(0,l, .928995))+.5(0))] 
+.360935(.3)[.5(1 - bin(0,2,.897144))+.5(l - bin(l,2,.897144))] 
+.548509(.3)[.5(1 - bin(0,3,.928995))+.5(l - bin(l,3, .928995))] =.27254 

Ps{20 min) =Ps{G) + /?X5) + A(10) + M15) = .84837. 

The expected number of missions executed is .84837. 
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APPENDIX E. COMPUTATIONS FOR THE MARKOV MODEL FOR 

CRUSADER AVAILABILITY 


This appendix provides the computations for determining, oi{H{t), t), the long-run 
availability of a Crusader pair. The value is used in the MMC. A numerical example of 
the computations is included, also. 

A. COMPUTATIONS 

The long-run proportion of time a Crusader pair is available given Crusader 
pairs are alive is; 


= limP(j)”^ (E.l) 


where P{y) is defined in (4.1). 
Let 


n->oo ■> 


(E.2) 


be the limiting distribution for the Markov chain with transition matrix (4.1). Therefore, 
|;r(y);y = satisfies the matrix equations 

nP{y) = Tt 

S^0)=1- <E.3),(E.4) 

j 

where n is a row vector. To solve the equations numerically rewrite (E.3) as follows; 

0 = ;r(/-P) (E.5) 

where / is an identity matrix with the same dimensions as P and 0 is a row vector of 


59 



zeros. Replace the third column of (/ - p) with 1 ’s. Solve with matrix operations for n 


as follows: 




Pn.E'-./P [a +a] 


c 

-Pba 

-Pra 


-PmlLrJ 


c=\ 




\yJ 


Pb 


Pba 

0 


(E.6) 


Solve for n{A),n{B),7r(R). The long run availability ;r(>’) = 7 r{A). 


It should be recognized that if the parameters in the model are time-dependent then 
there will not be a long-run availability. In this case, solve for 7 i(y) by matrix 
multiplication of (4.1) and use the value in (3.1). 

B. NUMERICAL EXAMPLE 


A detailed definition of the parameters which follow is included in Appendix B. 
L Parameters 


Pmit) 

■3; 

ri 

.5; 

ri 

.5; 

y 

= 1,2,3 pairs; 

Pb 

= .1111 computed in Chapter V; 

Pba 

= .347222 computed in Chapter V; 

Pr 

= .00245 computed in Chapter V; 

Pra 

= .01292 computed in Chapter V. 


2. The Transition Probability Matrix 

For the case where the number of Crusader pairs alive is 2: 


P(2) = 


1-.3(.5(1 / 2)(.l 111+.00245)+.5(2 / 2)(.l 111+.00245)) 
.347222 
.01292 


.3(.5(l/2).llll+.5(2/2).Ull) 

.65278 

0 


.3(.5(1 / 2).00245+.5(2 / 2).00245) 
0 

.98708 
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.02555 -0250 1 

/-P(2)= -.347222 .65278 ] 

-.01292 0 1 

[o 0 \\l-m]'=[ k(A,) ;r{J},)]; 

[o 0 l|/-?(2)]^' =[,897144 .064588 . 038268] 

P(l) and P{3) are computed similarly to provide the long-run probabilities for Crusader 
pairs available. The cumulative computations result in the following matrix: 

Avail Resupply Repair 
.928995 .044588 .026418 
.897144 .064588 .038268 
.928995 .044588 .026418 

where a{H{t) = 2; t) = .897144 in the MMC. 



1 

Crusaders Alive 2 
3 
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APPENDIX F. INSTRUCTIONAL DOCUMENTATION FOR THE CBSI 


This appendix presents the instructions included in a “readme.doc” file for the 

CBSI. 

1. Enter a battletime into the input dialog box in whole, integer hours from 1 to 96. 

2. Enter all input in appropriate units as stated on the CBSI input page. 

3. If an error occurs as a result of improper input, select “End” from the “pop up” 
dialog box then check to ensure all input is entered correctly. 

4. The probabilities of a type c mission arriving must sum to 1. 

5. The sheet with the name “GraphicsPage” is a blank sheet which may be used to 
construct graphs. 

6. Enter all input data, as desired, then use the “RECOMPUTE” or “FINITE-TIME” 
button to invoke calculations. 

7. Each sheet in the CBSI is named for organizational purposes. The name identifies 
the component of the CBSI which appears on the specified sheet. The sheet names 
are as follows: 

a. Input_Output: user input and output page 

b. GraphicsPage: page available for graphics construction 

c. IntComp: intermediate computations page 

d. FiniteTimeComp: finite-time computations page 

e. VBFiniteTime: Visual Basic program module 

8. The “VBFiniteTime” sheet is protected to prevent the program firom being 
changed, inadvertently. The password to change the protection status is: “cbsi” 
(no quotations). The procedure to “unprotect” the sheets is as follows: 

a. select the “VBFiniteTime” sheet; 

b. select the Tools menu; 

c. choose Protection; 

d. select Unprotect Worksheet. 
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APPENDIX G. CRUSADER BATTALION SPREADSHEET IMPLEMENTATION 

COMPUTATIONAL PAGE EXAMPLE 


This appendix presents an example of an intermediate computational page from the CBSI for 
the “base case” Crusader data file. Figure 12 displays the results of intermediate calculations for the 
supporting models of the MMC. These results are used in the final calculations for the MMC. 


Intermediate Computational Results for the CBSI 


Time for Mission = 


0.4 (nf/theta) 


Probability 1 Crusader Pair 
Engaged in Fire Mission of 
Target Type c is Killed: 

Pair c =1 

Platoon c =2 

c =3 

Battery c =4 

c=5 
c=6 
c=7 

Battery c -8 

c=9 
c =10 
c=11 

Battalion c =12 


0.112639 

0.199903 

0.266079 

0.314811 

0.349189 

0.371828 

0.384936 

0.390374 

0.389703 

0.384231 

0.375046 

0.363056 


Probability Crusader Pair 
Needs Resupply = 

Probability Crusader Pair 
Completes Resupply = 

Probability Crusader Pair 
Needs Repair = 

Probability Crusader Pair 
Completes Repair = 


hk(c) 


0.112639 


jProbability 


Red Acquires Blue = 


0.91 


Shoot and Scoot Time = 


Rate at Which Red Engages 
Blue = 


1.906325 


0.069315 


0.111111 


0.347222 


0.002451 


0.01292 


Figure 12. Intermediate Computational Results from the CBSI 
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APPENDIX H. VISUAL BASIC FINITE-TIME COMPUTATIONS CODE 


Sub FiniteTimeO 

Dim time As String 
Dim counter As Integer 
i = 4 
j = 16 

sumexecute = 0 
sumlost = 0 

time = InputBoxC'Provide number of hours of battle (min 1: max 96)") 

steps = Clnt(time) * 12 

Sheets("Input_Output"). Select 

Cells(88, 4).Select 

Selection. Value = time 

Calculate 

'Computes the P(2) matrix 
Sheets("FiniteTimeComp"). Select 
Range("a4:ml6").Select 

Selection.FormulaArray = "=mmult(transmatrix,transmatrix)" 

'Computes the P(3) to P(n) matrix and cummulates expected missions 
'lost and executed 
For counter = 1 To steps - 3 
i = i+ 14 
j = j + 14 

Range(Cells(i, 1), Cells(j, 13)).Select 

Selection.FormulaArray = "=mmult(r[-14]c:r[-2]c[12],transmatrix)" 

Cells(i + 12, 14). Select 

Selection.FormulaArray = "=mmult(rc[-13 ] ;rc[-l ],probmsnlostmatrix)" 
sumlost = sumlost + Selection. Value 

Cells(i + 12, 15). Select 

Selection.FormulaArray = "=mmult(rc[-14]:rc[-2],probmsnexecutematrix)" 
sumexecute = sumexecute + Selection.Value 

Next counter 

'Computes partial missions lost and executed in battle time 
'(does not include lost/executed at t=0,5,10 minutes) 

Cells(39, 17).Select 
Selection. Value = sumlost 
Cells(41, 17).Select 
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Selection. Value = sumexecute 


'Computes expected value and variance of Crusaders alive 
Range(Cells(i + 12, 1), Cells(j, 13)).Select 
Names.Add Name:="temp", RefersToRlCl;=Selection 
Cells(24, 15). Select 

Selection.FormulaArray = "=mmult(temp,Numbermatrix)" 
Cells(19, 15). Select 

Selection.FormulaArray = "=mmult(temp,Squarematrix)" 
Cells(22, 15). Select 
Selection.Formula = "=ol9-o24^2" 

'Recomputes spreadsheet values 
Calculate 

Sheets("Input_Output"). Select 
Cells(96, 6).Select 

End Sub 
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APPENDIX I. /MLAZ)/NCBSI INPUT PAGE 


Crusader Battalion Spreadsheet Implementation 

* 12 Crusader Pairs ( 24 Total Crusaders ) 

* Pair to Battalion Fire Missions 


User Provided Input: 


Probability Fire Mission Arrives: 



0.3 

Median Time until Red Engages Blue = 

10 

Probability Fire Mission is Target ( min ) 


Type c Requiring c 

Crusader 



Pairs to Fire: 


Number of Rounds Fired per Mission = 

4 

Pair Mission 

0.01 



Platoon Mission 

.©.75 

Rate of Fire (rounds/min) = 

: . 4 





Battery Mission 

•••0,1 

Number of Missions Fired Before 




Survivability Move = 

1 


om 




wz 

Resupply Ammunition at Percentage 


Battery Mission 

0.002 

Expended (decimal percent) = 

0.6 


0.002 




0.002 

Ammunition Transfer Rate (rnds/min)= 

1.714 


• 0,002 



Battalion Mission 

•.:0,1 

Mean Time Between Essential 




Function Failures (hours) = 

■•^^.21 

Scoot Speed(m/s)= 

^ -0,22 



(750m/90s»8.3m/s) 


Mean Time to Repair plus ALDT = 

6.45 


• * > 

(hours) 


Counterbattery 


Probability of Acquisition = 

.0.7 

Footprint(meters)= 

.080 

of Crusaders by enemy 




Sensors as a Result of 


Crusader Payload = 

••■ 39 

Firing 
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